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Abstract 

Ballistic annihilation kinetics for a multi-velocity one-dimensional ideal gas 
is studied in the framework of an exact analytic approach. For an initial 
symmetric three-velocity distribution, the problem can be solved exactly and 
it is shown that different regimes exist depending on the initial fraction of 
particles at rest. Extension to the case of a n- velocity distribution is discussed. 

PACS numbers: 82-20.-W, 05.20Dd, 05.40-Fj 



I. INTRODUCTION. 



Ballistically-controUed reactions provide simple examples of non-equilibrium systems 
with complex kinetics. In the one-dimensional case, one considers point particles which 
move freely, with a given velocity. When two particles collide they instantaneously annihi- 
late each other and disappear from the system. The system with only two possible velocities 
+c or — c has been studied in a pioneering work by Elskens and Frisch . Using combinato- 
rial analysis, they showed that the density of particles was decreasing according to a power 
law (t^^/^) in the case of a symmetric initial velocity distribution. Later, Krug and Spohn P] 
obtained independently similar results. Recently, Redner and co-workers have studied 
(essentially numerically) the problem with more general velocity distributions. However, no 
rigorous analytical approach has been elaborated until one of us (JP) developed a formalism 
in which the annihilation dynamics reduced exactly to a single closed equation for the two- 
particle conditional probability 0. The purpose of this work is to present a method which 
permits to solve this evolution equation for discrete velocity distributions. More precisely, 
we shall consider the three- velocity case where the initial velocity distribution 0(f ; t = 0) is 
given by 

(j){v] 0) = p+5{v - c) + Pq5{v) + p^5{v + c) (1) 

with J9+ = (symmetric case), and p+ + po + p_ = 1. 

For the ballistic motion, the collision frequency between two particles is proportional to 
their relative velocity. Then, collisions involving particles at rest occur less often than the 
ones involving moving particles. Indeed, the collisions between one particle of velocity +c 
and one of velocity — c (relative velocity 2c) occur twice more often than those between a 
particle at rest and a moving one (relative velocity c). As a consequence, one can expect that 
in one half of the total number of collisions participate two moving particles and the other 
half involves one particle at rest. Thus, in mean, only one quarter of particles annihilated 
in collisions are at rest and the three remaining quarters are formed by moving particles. 
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This simple heuristic argument suggests that the values Pq = j, pX = P- = ^ must play 
a special role: if po is less than p^, the system will asymptotically behave as in the two- 
velocity case (the stationary particles will disappear before the annihilation of the moving 
ones). On the contrary, if po is greater than Pq, the moving particles will disappear first, and 
asymptotically only particles at rest will be left. In the limiting case po = p^ the moving 
particles and the ones at rest disappear at the same rate. It will turn out that this intuitive 
argument is exact. We shall show it using the approach presented in [^], which permits to 
solve this three-velocity model exactly. 

This paper is organized as follows. In section |I|, we introduce the model and recall 
the main steps of the rigorous method [^. An important quantity is S{v;t), the survival 
probability until the time t of a particle moving with velocity v. In section pT| , we study the 
behavior of S{v; t) and compute the density and the time dependent velocity distribution in 
the asymptotic regime t — > oo. The value = i appears to be a critical point, separating 
differents kinetic behaviors. In addition, it is shown that the dynamics is incompatible 
with a Boltzmann-like approximation due to the appearance of strong velocity correlations 
between particles which are nearest neighbors. The case of a general discrete multi-velocity 
distribution is examined in section Finally, concluding remarks are given in section [V|. 

II. THE THREE- VELOCITY MODEL. 

We assume that initially the particles are uniformly distributed in space, according to the 
Poisson law, without any correlations between their velocities. Note that other distributions 
than Poisson could be considered as long as one is dealing with a renewal process. The 
process is thus translationally invariant and each particle has initially the same probability 
density 0(f ; 0) (given by equation (|1|)) to move with velocity v. 

The important characteristic of the annihilation dynamics is that only those particles 
which suffered no collisions are present in the system. They are thus found on their free 
trajectories. Accordingly, the key quantity of the theory is the survival probability S{v; t) 
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which has the product structure 

S{v;t) = S''{v;t)S''{-v;t) , (2) 

where S^{v; t) is the probabihty for the absence of colhsion with the right neighbor. As 
S^{-c; t) = 1, for t > 0, we find 

5(+c;t) = 5(-c;t) = 5^(+c;t) , (3) 

Initially, S^{+c; 0) = S^{0; 0) = 1. 

The density of particles with velocity v at time t is given by: 

a{v;t) = aS{v;t)(f){v;0) , (5) 

where a is the initial density. 

Another important quantity in the approach |^ is the distribution of nearest neighbors. 
Suppose that at time t there is a particle moving with velocity v. We denote by mIO, v; t), 
the conditional probability density for finding its nearest neighbor to the right at a distance 
X > 0, with velocity u (for non-homogeneous systems the probability may depend upon the 
positions of both particles, rather than on their relative distance only). As we are dealing 
here with symmetric initial velocity distribution, the following relation holds for /i: 

; 0) S{v; t) u\0, v; t) = (f>(—u] 0) S{—u; t) — f |0, — m; t) . (6) 

The initial condition is given by 

m|0, t = 0) = d{x) ae-'"'(f){u; 0) , (7) 

where 9{x) is the usual Heaviside unit step function. 
We also define the density at contact by 

fj.{0~^ ,u\0,v]t) = lim fi{x,u\0,v]t) . (8) 

a;— >0 
x>0 
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5(0; t) = 5^(0; t) 



This quantity plays a particular role, because it determines the density of precollisional 
configurations. 

The time evolution of S^{v; t) is given by: 

S\+c; t) = exp |- j^^dr [c/i(0+, 0|0, +c; r) + 2c/i(0+, -c|0, +c; r)] | (9) 
S^{<d]t) = exp|-^Wc/x(0+,-c|0,0;r)| (10) 



and the two-particle conditional probability density obeys to the following closed equation 
\dt^ '''''' ^ S^{v^,t) S^^iv^-^t) 



+ ^^21 + 7757 TTT - -TT > IJ'{z,V2\Vi;t) - V21 H{0~^ , V2\0, Vi, t) 



dVsJ dVifl{z,V3\Vi;t) fl{z,V2\v4;t)vuO{v34) f^{0^,V4\0,V3]t) , (11) 

where the dot stands for the time derivative, V21 = V2 — vi, and we use the definition 

rex 

jl{z,u\v;t) = J ^dxe~^^fi{x,u\0,v;t) . (12) 

It should be stressed that this remarkable property of getting a rigorous closed equation 
for the distribution /i comes from the main characteristic of the annihilation process: only 
those particles which have not interacted could survive until time t. For more details see 
reference P] . We can now start to solve equation (|l^) for our symmetric three velocity case. 
Equation ([TlD can be rewritten in a matrix form: 



J\f{z; t) + C{z; t) = Af{z; t) C{z; t) Af{z; t) (13) 
where A'' and C are two 3x3 matrix: 

Ar(^;t) = (%(^;t)).^.^^2,3 (14) 

with 
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and 



with 



N2i{z 

N23{Z 

N32{Z 
N33iz 



t) = fi{z, +c| + c; t) 

t) = ji{z,0\+c;t)e-'^'S'^{+c;t)/S'^{0;t) 

t) = jl{z,-c\+c;t)e-^^'*S^{+c;t) 

t) = fl{z,+c\0;t)e"^S'^{0;t)/S'^{+c;t) 

t) = jl{z,0\0;t) 

t) = jl{z,-c\Q;t)e-^^ S^{Q;t) 

t) = jl{z,+c\-c;t)e'^''''/S^{+c;t) 

t) = jl{z,0\-c;t)e'''yS^{0;t) 

t) — jl{z, —c\ — c; t) 



^0 Cu{z-t) Cn{z-t)^ 



C{z;t) 



C23{z;t) 




Cu{z;t) = cfi{0+,0\0,+c;t)e-"''S^{+c;t)/S^{0;t) 

Ci3{z;t) = 2c//(0+,-c|0,+c;i)e-2^=*5^(+c;i) 

C23{z;t) = c/i(0+,-c|0,0;t)e-^^*>5«(0;t) . 
The initial values of and C are: 



J\f{z; 0) 



z + a 



I \ 

P+ Po P+ 

P+ Po P+ 



\P+ Po P+ j 
^ ccrpo ^cap-^ ^ 



C{z;0) 



C(Tp+ 
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In spite of the fact that equation ([T3D looks very simple, it should be noticed that it is a 
matrix differential equation for which the general solution is not known. However, the par- 
ticular structure of C makes it solvable. The method of solution is described in appendix 
We find 



iVii 
Nu 

N23 

N32 
N33 



{z;t) 

{z;t) 
(z-t) 
{z;t) 
{z;t) 
(z-t) 
(z-t) 
(z-t) 



D{z-t) 

_Poo_ 
D{z-t) 

D{z-t) 

D{z-t) 

_PoO_ 
D{z-t) 

D{z-t) 

D{z-t) 

D{z-t) 

D{z-t) 



1 



_Po_ 

n_ 

2p+ 



dTC2?,{z\ t) 



dTC2?,{z\ t) 



p+ 



[ dTC2'i{z]T) 

Jo 

l^drC23iz;r)j') - J^drC,siz;r) 



(20) 



Po 
2p+ 



dTC23{z;T) 



dTC23(z: T 



Po 
2p+ 



dTC2,3{z; t) 



where 



D{z;t) = z + a-2poa f dr C23{z] t) - p+a f dTCri{z]T) 

Jo Jo 



(21) 



The Laplace transformed conditional probabilities fl are obtained from equations (jTH]). 
Thus we found the exact solution to the kinetic equation (0). However, this is as yet an 
implicit solution, that is to say, the values of /2 are expressed in terms of the density at 
contact fi{0^ ,u\0,v]t) {u < v). To obtain some information about physical quantities, we 
have to solve the consistency equations which express the three non-vanishing densities at 
contact in terms of function jl: 



7+ioo 



dz 
2i7r 



fi{z, u\v] t) {u < v). 



(22) 



7— JOO 



Up to this point, the genrahzation to a non-symmetric three- velocity distribution is 
straightforward. 



III. LONG TIME BEHAVIOR OF THE SYSTEM. 

Instead of writing the consistency equation for the density at contact, we shall write 
them for the survival probabilities. Remembering that only two densities at contact are 
independent (see equation (^)) and that (from equations and (pUD) 

S^(0; t) = -c/i(0, -c|0, 0; t) S^{0; t) (23) 
S^{+c; t) = [-c/i(0, 0|0, +c; t) - 2cn{0, -c|0, +c; t)] S^{+c; t) , (24) 



we obtain, using equation (^4 ), the following two equations: 

^^(0;t) = -cp+cr / ^e^^* ' ^ (25) 



7—100 



2i7r D{z;t) 



5^(+c; t) = ^S\0- 1)5«(0; t) (26) 
P+ 

JqcIt C23iz; t) 



7+ioo 1 _ PO 

-2cp+a / —e 



2m D{z;t) 

7 — JCXD 

Solving this two equations seems to be extremely difficult, principally because of the time 
dependence in D{z] t) and in the integrals over r. However it can be shown (see appendix jB]) 
that the only relevant time dependence occurs in the exponential factor. By this we mean 
that equations (|25|) and ( ^6]) can be exactly replaced by the following two relations: 

= -cp^a / ^e-* Z . r (27) 



2z7r D{z; 00) 

7— iOO 

S^(+c; t) = —5^(0; t)S^(0; t) (28) 
P+ 

^+;- ^ fl - U^dTC23(z;r)?Y 



2m D{z; 00) 

7—100 

This remarkable fact allows us to use again the Laplace transform to suppress the inverse 
Laplace integral and arrive at much simpler relations: 
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T{z) 



1 - i;T\z) 



z + a + 2poaT{z) + p+alJ (z) 



(29) 



U{z) = dte 
Jo 



-2zct 



5^(+c;t)-— 5^(0;t)5^(0;t) 
l-^T\z'" 



z + a + 2pQaT{z) + p+aU{z) 
where equations (|17[) and ( pT]) have been used. From this system of equations it follows that 
T{z) satisfies the quartic equation 



Po T\z) - {2po + 1) T^{z) - 2 (^^ + l) T{z) + po - 1 = , 



U{z) being given by 



U{z) = T{z) 



Po rr.3 



2p, 



(31) 



(32) 



It is hopeless to try to Laplace invert the solution of equation (^) to find S^{0] t). However, 
the asymptotic regime of the survival probability can be determined by considering T{z) for 
z in the neighborhood of 0. For z = 0, one knows that T(0) is given by 

/•OO 

r(0)=/ dtS^{0;t) = S^iO,oo)-l (33) 
Jo 

and that 5^(0,00) G [0, 1]. This condition allows us to identify the physically acceptable 
solution of the above quartic equation (^). We have to distinguish between two different 
situations: 



Po < 



The only acceptable solution is T(0) = —1, corresponding to an asymptotically empty 
stationary state S^{+C] oo) = 5*^(0; oo) = 0. 



Po>\ 

We have two possible solutions: T(0) = —1 and T(0) 



However, S^{0] t) 



is a continuous decreasing function of time with initial value 1. Thus the stationary 
value 2 — will be reached first. Accordingly, T(0) = 1 ^ is the relevant value. 
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We can now study the behavior of T{z) near zero, by introducing the quantity e{z) = 
T{z) — T(0) in the quartic equation (^). Five cases have to be distinguished: 

• Po = 

This is the bimodal velocity distribution aheady investigated "^^^ quartic equa- 
tion ( piD simphfies to the following quadratic equation for U{z): 

U\z)+2{- + l)U{z) + l = . (34) 
a 

By solving and Laplace inverting it, one recovers the Elskens and Frisch results [|1|. 

• < po < i 

Now e{z) = T{z) + 1 obeys the quartic equation: 

Poe\z) - 4poe^{z) + (4po - l)e'(^) - 2-{e{z) - 1) = . (35) 

a 

As we are interested in the limit z ^ 0, and as lim^^oe(-2) = 0, we can neglect the 
terms of order e"^ and e'^ with respect to e^, and e with respect to 1. By Laplace 
inverting and integrating, we obtain the following asymptotic behavior {t oo): 



For t oo, the densities are then given by: 

"'"■"^ (l-4po)art {^ + °(l^-^P°l'>'l")} • 

= \/(i-Po)£j{l+0(|l -4Po|-'"M-"')} ■ (39) 

The corrections to scaling have been obtained by a careful study of the solution of the 



quartic equation (^). 

The prefactor {Apo — 1) in the term in equation ( |55D indicates that for po near to | 
we should expect important crossover effect. This is confirmed by the amplitude of the 
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correction to scaling for a{0;t) which contains a term proportional to (1 — 4po)~^- As 
a consequence, it may be very difficult to extract the true asymptotic behavior from 
experimental or numerical data. However, we see that for very long times the system 
is driven towards the bimodal case. 



Po = 4 

Here again, e{z) = T{z) + 1. The coefficient of the term vanishes and the term of 
order cannot be neglected. We are thus left with the cubic equation: 

e^(^) - 2- = 0, when ^ ^ (40) 

(7 

whose physically relevant solution leads to the asymptotic behavior (for t — > oo): 



r(2/3) \catj 
1 / 2^3 Y 



3 vr(2/3) ) r(i/3) 

For the density one gets {t oo): 

/ 2 



\catj 



a{0;t) 



a 



412(2/3) \cafj 



3a 

a{+c]t) = — 



1 / 21/3 y 1 



+ 



3 vr(2/3) J r(i/3) 



i-f 



{l + O {[cat]-^')} . 



In addition. 



^ _ 1 ^ 3 r^(2/3) ^ 

a{0;t) 2 2 22/3r(l/3) 



(41) 
(42) 



(43) 
(44) 



(45) 



Now e{z) — T{z) — 1 + ^=- The terms of order and can be neglected with respect 



to the term and variable z with respect to 1. Hence: 



e{z) 



Po(5 - 2^) 



(46) 



2(1 - v^) 



\po(2VP^-l)(5-2vp^) 



[2VP^-1)= 



z 

+ - 



2(1- VW)(5-2VP^) a 



1/2 
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which leads to an exponential asymptotic behavior for S^{0; t) —2 + and also for 
S^{+c;t). Thus, for i ^ oo 

(7(0; t)^a (2^ - if + 2a Aj—^ , (47) 

Q—caut 

-(+^;^)--^(^' (48) 

where 



A^ 



2po(l- V^)(5-2^) 



u 



\ 71(2^0-1) (2^^-1)2' 

(2^-1)3 
2po(l- V^)(5-2v^) • 



In addition, 



hm "<''^'>-"'°^°°) = 2 . (49) 

indicating that in the long time regime the collisions between pairs of moving particles 
disappear, and only collisions involving a particle at rest and a moving one can be 
observed. 

• Po = 1 

In this trivial limiting case, in which there is no collisions, the physical solution is 
T{z) = 0, and the density remains constant. 

As anticipated on heuristic arguments, we have shown that po = | is playing a special 
role separating the different kinetic regimes. 

Another important quantity that can be computed exactly concerns the correlations be- 
tween the velocities of the colliding particles. Although at time t = these velocities are 
uncorrelated, annihilation dynamics creates strong correlations between them during the 
time evolution. This has the important consequence to exclude a Boltzmann-like approx- 
imation. It is clearly seen on w{y\t\ the mean velocity of the right nearest neighbor of a 
particle moving with velocity v. 
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w{v; t) = c/i(0, +c\v; t) — c/i(0, — c|f ; t) . 



(50) 



Using equations ( poD we obtain for example: 



w{—c; t) 



D{0;t) 



with D{0;t) given by equation ( pTD and 6*23(0; t) by equation (0). Our results show that 
pairs of nearest neighbors have the tendency to align their velocities and propagate in the 
same direction. Indeed, in the limit t ^ 00, one finds for Po < j- 



for Po 



w{v; 00) 



±c 
c(V2 



1 -4po 
l-2po 



if f = ±c 
if = 0, 



(51) 



w f ; 00 = V 



(52) 



and for Po > i 



w[v:oo\ 



-1/2 



(53) 



IV. GENERAL DISCRETE VELOCITY DISTRIBUTION 

In this section, we want to consider the case of a general discrete velocity distribution. 
Suppose that the initial velocity distribution is given by 

n 

(p{v;0) = Y,Pk5{v -Vk) , (54) 

k=l 

with 

n 
k=l 

As for the three-velocity case, we assume initially random spatial distribution and no corre- 
lations between the velocities. We can generalize some of our formulae in a straightforward 
way: the probability for the absence of collisions with the right neighbor is given by: 
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S'^iv; t) = exp y^rj^iv - Vk) e{v - t;fc)/x(0+, Vk\0, v;t)^ . (55) 
If we define = 1, . . . , n): 

Nij{z;t) = fM{z,Vj\vi;t)expi^-zVijt + J^dTY^^[Bik{T) -5jfc(r)]| (56) 
C,,{z- 1) = B,,{t) exp I -zv,,t + f^dr J2 [Bik{r) - 5,fc(r)] | (57) 

where 

Bij{t) = Vij9{vij) n{0,Vj\0,Vi;t) , 

the evolution equation for the Laplace transformed conditional probability fl (equation |TT]) 
can be put in the same form as equation ([13|), but with a n x n matrix. Once again, the 
particular structure of the matrix C will help us in solving this equation: assuming that 
A/'(-2; t) — X is invertible (this can be explicitly checked for t = and verified a posteriori for 
later times), where X is the n x n unit matrix, we define 

V{z;t) = mz;t)-I)-' + ^ . (58) 

Then, equation (|T3|) reads 

^(z; t) = -Viz; t)C{z; t) - C{z; t)V{z; t) . (59) 

One way to solve this equation is to remember that Cij{z; t) = for i > j {i, j = 1, . . . ,n). 
Hence, equation (|59|) has the following structure: 

j-l n 

- {z; t) = J2 Pik Ckj + E Cik Pkj (60) 

k=l k=i+l 

and thus when i = n, j = 1 it takes the form 

P„,i(2;;t) = . 
In addition, the [n — 1,1) and {n, 2) equations are 
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— Pn~l,l(yZ; t) — Cn-l,n{z'i t)Pn,l{z] t) 

-Pn,2{z;t) = Ci,2{z;t)Pn,i{z;t) , 

and can be solved once Pn,i{z; t) is known. One can easily see that the equations for {n—2, 1), 
{n — 1, 2) and (n, 3) are expressed in term of Pn,i, Pn-i,i and P„,2, allowing us to complete the 
solution. The process is iterated until solving equation and determining entirely the 

matrix V{z;t). The Laplace transformed conditional probabilities are obtained straight- 
forward by inverting the formula (^8]). Again we have an implicit solution and we have to 
solve the n{n — l)/2 consistency equations (see equation (|2^)) to obtain physical information 
on the system. Although detailed calculations can be very tedious, there should not be any 
conceptual difficulties. 

V. CONCLUDING REMARKS 

We have shown in this paper that new exact predictions on ballistic annihilation kinetics 
can be obtained in the framework of the new approach [Q. It is remarkable that, for discrete 
velocity distributions, the non-linear integro-differential equation governing the two particle 
conditional probability can be solved exactly. 

For po < ^, not only the leading power laws are obtained in the asymptotic regime, 
t — > oo, but also the amplitudes and the corrections to scaling. In particular, we have shown 
that important corrections to scaling occur when po is in the neighborhood of |. For Po > \, 
exponential behavior is obtained for long times. The results of our numerical simulations 
and those of Redner et al. are well explained by our exact theory. In addition, it is 
clearly shown that Boltzmann-like approximations fail, because the annihilation dynamics 
favorizes configurations in which the nearest neighbors have the same velocity. 

Several extensions of this work are possible. In particular, the case of continuous velocity 
distributions for which qualitatively different behavior may be expected is under investiga- 
tion. 
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APPENDIX A: ON THE MATRIX EQUATION 



The solution of the matrix equation (|T3D is presented here. To solve it, we can write M 



as 



N{z-t) 



(Al) 



where Mi is a 1 x 2 matrix, A/2 a 1 x 1 , A/3 a 2 x 2 and Ma a 2 x 1 . In this decomposition, 
C is given by 



C{z-t) 



^ Si{z;t) S2{z;t)^ 
y £3iz;t) ^ 



(A2) 



where 



^iiz;t) = (0 Ci2{z;t) ) , 
£2{z;t) = ^Cis{z;t)) , 
( C,..(z;t) ^ 

S3{z-t) - 



V 



C23{z;t) 


However, to perform the matrix product MCM, we have to write C in a different way, 
namely: 



C{z;t) 



( £{z-t)^ 







(A3) 



where £^ is a 2 x 2 matrix 



£{z-t) 



^ Ci,2{z-t) Ci,3{z;t)^ 



\ 



C2,3iz;t) J 

With this decomposition, the matrix equation (^) becomes 
/ . 

Mi{z;t) M2{z;t) 
M3{z;t) M^iz-t) 



I 



\ 
i 



E^^z-t) 



(A4) 



Mx{z; t)S{z; t)M3{z; t) Mi{z] t)S{z; t)Mi{z] t) 



\ 



liz;t)S{z;t)M3iz;t) M3iz;t)S{z-t)Mi{z;t) j 
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Although TV's is not invertible, we can find a solution for the equation 

U^{z; t) = Mi{z; t)£{z; t)Afs{z; t) (A5) 

using the following Ansatz: 

Afsiz; t) = A{z) [I - B{z- t)A{z)]-^ (A6) 

where X is the unit 2x2 matrix. Once we have obtained 7V3, we can readily solve the three 
other matrix equations and thus find M . 

APPENDIX B: ON THE CONSISTENCY EQUATIONS 

In this appendix we show how to justify equation (|27|) and ( p8D starting from equa- 
tion (0) and (|6|). 

A way of demonstration could be to simply subtract equations ( pTf ) and (|25| ) and to 



verify by explicit integration in the complex plane, that the result is zero. Instead of going 
through all this tedious algebra, we prefer to give a simple physical argument. 
We define Ai{z;t) by 

+c| + c; t) = p+aAi{z; t) . (Bl) 



From the kinetic equation ([TTD , we have 



3 3 

jj,{z,+c\+ c]t) = ^^jx{z,Vi\+c]t)jx{z,+c\vft)vijd{vij)^,{Q'^,Vj\Q,Vi]t) . (B2) 
i=i j=i 

Now, the right hand side of equation ( p2D represents the variation of jl due to the mu- 
tual annihilation of particles originally separating the pair (+c, +c) (see reference for 
the interpretation of different terms in equation (p!T|)). It follows that the inverse Laplace 
transform 

7+ioo 

Ai(x;t)= / iie^Mi(z;t) (B3) 

J 2t7T 
7— ioo 
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represents the probability that the particles present initially in an interval of length x, 
separating two particles with velocity +c, disappear through ballistic annihilation before 
time t (the prefactor p+a in equation (|B1| ) is the initial density of the right neighbors 
moving with velocity +c). As the colliding pairs move at least with a relative velocity c, all 
sequences of encounters contributing to Ai{x;t) are accomplished at the moment t* = x/c. 
Therefore Ai{x; t) does not depend on time (for fixed x) in the region t > t*. 
Now we remember that 



1 - 



PO 
2p 



JodTC23iz;T 



2 



Mz;t) = — , (B4) 



D{z;t) 

hence the equation (^5|) reads: 



7+joo 

S^{0;t) = -cp+a I -±e^-'A,{z-t) . (B5) 

7— joo 

In the term on the right hand side of equation (BH), the inverse Laplace transformation 



yields the value of function Ai{x;t) at the point x = ct. So we can replace Ai{z;t) in 
equation (p5| ) by Ai{z;oo), without changing the value of the integral. Hence we obtain, 
equation (|27|) . 

The equation (^6|) can be handled in the same manner: we first define 



P+(tA2{z; t) = fi{z, -c\ + c; t) e"'""* S^{+c; t) + f dt Cn{z; t) . (B6) 

Jo 

Then, by inspecting the evolution equation for /i(z, — c| + c; t) we remark that ^2(2;; t), the 
inverse Laplace transform of A2{z]t), represents the probability that all particles present 
initially in an interval of length x, separating one particle of velocity +c on the left and one 
of velocity — c on the right, disappear through ballistic annihilation before time t. It follows 
that A2{x; t) does not depend on the time (for fixed x) in the region t > t* /2 (as the relative 
velocity between the colliding pairs must be 2c). Finally, the expression 



1 _ PO 

2p4 



JodT C23iz; t) 



leads to the expecting conclusion. 
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